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北斗 三 号 卫星 精密 定 轨 中 的 光 压 模型 研究 : 


SA KHL WEE REF 


(河南 理工 大 学 测绘 与 国土 信息 工程 学 院 焦作 454000) 


摘要 对 于 在 轨 运 行 的 BDS (BeiDou Navigation Satellite System) 卫 星 , 太阳 光 压 是 作 | 


在 卫星 上 主要 的 非 


引力 摄 动 . 受 多 种 因素 的 影响 , 太阳 光 压 摄 动力 难以 精确 建 模 , 是 BDS 卫 星 精密 定 轨 和 轨道 预报 过 程 中 重要 的 误 


差 来 源 . 由 于 ECOMC (Empirical CODE Orbit Model 1 and 2 Combined) Z2 3i Bii T E 


COMI (Empirical 


CODE Orbit Model 1) 和 ECOM2 (Empirical CODE Orbit Model 2) 模 型 的 特点 , 在 模型 中 引入 了 较 多 的 待 


© 对 该 模型 的 光 压 参数 进行 时 间 序列 分 析 ， 


~ (HBR, 使 得 参数 之 间 存 在 强 相关 性 . 针对 ECOMC 模 型 的 这 一 缺陷 , 文中 收集 了 2019 年 1 月 
= 学 分 析 中 心 提供 的 BDS-3 卫 星 精 密 星 历 , 采用 动力 学 轨道 拟 合 方法 得 到 了 ECOMC 模 型 的 13 个 光 压 参数 . 通过 

分 别 给 出 了 BDS-3 IGSO (Inclined Geosynchronous Orbit) 和 MEO 
(Medium Earth Orbit) 卫 星光 压 模型 的 参数 选择 策略 . 并 利用 轨道 拟 合 和 轨道 预报 试验 , 验证 了 光 压 模型 参数 


至 2022 年 4 月 武汉 大 


时 , 也 能 够 提升 BDS-3 IGSO 和 MEO 卫 星 中 长 期 轨道 预报 的 精度 . 
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N 1 引言 方式 ， 以 改善 经 验光 压 模 型 的 定 轨 性 能 ; (ID 根据 
太 阳 光 压 摄 动 又 称 太阳 t J} JE (Solar Radia- BDS 卫 星 元 数据 或 者 由 可 校正 Box-Wing 模 型 估计 


选择 策略 的 合理 性 . 结果 表明 , 采用 改进 型 BCOMC 模 型 进行 BDS-3 IGSO 和 和 MEO 卫星 轨道 拟 合 的 效果 最 佳 , 同 


O 续 " 撞 击 "卫星 表面 所 产生 的 力 的 作用 . 对 于 在 轨 运 


= 行 的 BDS (BeiDou Navigation Satellite System) E 针对 (D 中 的 建 模 思 


"so E, 其 所 受 的 非 引 力 摄 动 主要 来 源 于 太阳 光 压 . 受 ” ”方法 , 提出 了 一 种 适 


tion Pressure，SRP) 摄 动 , 是 由 太阳 辐射 的 光子 连 。 的 卫星 几何 、 光 学 特性 参数 , 建立 BDS 卫 星 的 先 验 
光 压 模型 以 提升 光 压 模型 的 定 轨 性 能 

起 ,文献 图 采 用 实验 设计 的 

]TBDS GEO (Geostation- 


” ”卫星 姿态 以 及 卫星 表面 材料 特性 的 影响 , 太阳 光 压 ary Earth Orbpit) 卫 星 的 经 验光 压 模型 , 并 通过 地 基 


摄 动力 难以 精确 地 建 模 , 是 BDS 卫 星 精 密 定 轨 和 轨 GNSS (Global Navigation Satellite System) 技 术 定 


道 预报 重要 的 误差 源 [i-2. 轨 和 卫星 激光 测 距 (Satellite Laser Ranging, SLR) 


顾 BDS 卫 星光 压 摄 动 模型 的 建 模 历 程 , Hou 技术 定 轨 两 种 方式 , 验证 了 该 光 压 模型 的 性 能 . 文 
模 思 路 主要 分 为 : (在 经 验光 压 模型 的 基础 上 , 利 献身 通 过 分 析 ECOM1-9 (Empirical CODE Orbit 
dBDS 卫 星 的 精密 轨道 , 求解 经 验光 压 模型 的 参数 Model 1,， 9 表示 参数 数量 


) 模 型 的 光 压 参数 变化 规 


并 进行 分 析 , 确定 出 经 验光 压 模 型 参数 的 最 佳 组 合 t£. 给 出 了 BDS-2 各 类 轨道 卫星 的 ECOMI1 经 验 模 
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64 卷 天 


型 的 光 压 参数 选择 策略 ， 并 提出 了 在 ECOM1-5 模 
型 的 万 方向 ( 即 DYB 坐 标 系 的 刀 轴 方向 , 其 中 忆 轴 为 
星 至 太阳 方向 、Y 轴 为 卫星 至 地 球 方向 与 卫星 至 
太阳 方向 的 义 乘 、B 轴 为 D 与 Y 轴 的 又 乘 , B 轴 与 D 
和 Y 轴 构成 右手 系 ) 上 增加 周期 项 ， 以 吸收 BDS-2 


ChinaXiv 合 作 期 刊 


synchronous Orbit) 与 MEO 卫 星 定 轨 的 适用 性 . 

本 文 根 据 BDS-3 卫 星 的 元 数据 , 基于 Box-Wing 
模型 计算 了 BDS-3 卫 星 的 光 压 模型 先 验 力 . 采用 卫 
星 动力 学 轨道 拟 合 方法 得 到 ECOMC 模 型 的 13 个 光 
压 参 数 ， 并 对 BDS-3 IGSO 卫 星 和 MEO 卫 星 的 光 


卫星 零 偏 姿态 控制 模式 下 的 光 压 摄 动力 . 文献 [四 采 
几何 轨道 平滑 方法 反 演 了 ECOM1-9 模 型 的 光 压 
参数 , 并 对 参数 的 时 间 序 列 进行 了 分 析 . 通过 实测 
数据 进行 定 轨 验证 , 发 现 了 采用 只 解 算 5 个 光 压 参 
数 的 ECOM1 经 验光 压 模型 进行 BDS 卫 星 定 轨 的 精 
度 最 高 . 

文献 [6] 采 用 多 模 GNSS 试 验 (Multi-GNSS Ex- 
periment，MGEX) 测 站 的 BDS 观 测 数据 进行 定 轨 
Wik, 得 到 了 扩展 型 EECOM (Extended Empiri- 
cal CODE Orbit Model) 模 型 在 BDS-2 的 3 种 轨道 
类 型 卫星 定 轨 中 的 光 压 参数 选择 策略 . 在 GNSS 卫 
星 精 密 定 轨 中 , ECOM1 经 验 模 型 还 存在 着 未 被 模 
型 完全 吸收 的 光 压 摄 动力 , 尤其 是 对 于 非 立 方 体 
结构 的 GNSS 卫 星 斑 9. 因此 (IT) 中 的 建 模 思想 被 广 
泛 采 纳 , 考虑 到 目前 BDS-3 卫 星 的 元 数据 并 未 完全 
公开 , 文献 [10] 采 用 假设 的 BDS 卫 星 几 何 、 光 学 特 
性 参数 建立 了 具有 明确 物理 意义 的 先 验光 压 模 型 ， 
并 联合 ECOM1-5 模 型 , 进行 BDS 卫 星 精 密 定 轨 . 文 
献 [11-13] 根 据 可 校正 Box-Wing 模 型 估计 的 相关 参 
数 , 建立 了 先 验 Box-Wing 模 型 ， 并 结合 BCOM1 模 
型 进行 定 轨 实验 , 采用 该 模型 获得 了 比 纯 经 验光 压 
模型 更 好 的 定 轨 效 果 . 文献 14-15] 基 于 光 压 摄 动 的 
作用 机 理 、 卫 星 表面 的 光学 特性 参数 和 卫星 姿态 
控制 模式 等 , 建立 了 BDS 卫 星 综合 解析 模型 . 同时 ， 
根据 卫星 角 动 量 守恒 和 在 轨 遥 测 参数 的 变化 , 提出 
了 在 综合 解析 模型 基础 上 增加 常数 经 验 修正 项 的 
方法 , 使 事后 精密 定 轨 精度 达到 了 厘米 级 . 文献 [16] 
根据 BDS 卫 星 的 元 数据 ， 在 Box-Wing 模 型 的 基础 
上 建立 了 BDS-3 MEO (Medium Earth Orbit) LÆ 
的 先 验光 压 模型 , 并 通过 定 轨 验 证 和 SLR 技 术 定 轨 
检验 , 发 现 采 用 先 验 光 压 模型 可 以 显著 减 小 BDS-3 
MEO 卫 星 重 胎 弧 段 轨道 的 误差 和 SLR 定 轨 残 差 的 
标准 差 . 文献 [17] 从 轨道 内 、 外 符合 精度 等 方面 , 分 
it 了 ECOMI1-5 模 型 、 Box-Wing+TECOMI1-5 模 型 
和 拓展 版 ECOM 模 型 对 BDS IGSO (Inclined Geo- 


压 参数 进行 了 时 间 序 列 分 析 ， 分 别 给 出 了 适用 了 
BDS-3 IGSO 卫 星 和 MEO 卫 星 的 光 压 参数 选择 策 
Wk. 还 通过 轨道 拟 合 和 轨道 预报 试验 , 验证 了 BDS- 
3 IGSO 和 MEO 卫 星光 压 模型 参数 选择 策略 的 合理 
PE. 由 于 BDS-3 GEO 卫 星 精 密 星 历 没 有 足够 多 的 
FRA, 本 文 暂 不 将 BDS-3 GEO 卫 星 纳入 研究 范围 . 


2 BDS-3 卫 星光 压 模型 

为 满足 高 精度 定 轨 的 需要 ，Beutler 等 0 以 
ROCK 模型 bg 为 先 验 模型 结合 精密 轨道 数据 ， 建 
立 了 ECOM1-9 模 型 . 经 过 国内 外 大 量 学 者 的 研究 
改进 , 该 模型 出 现 了 众多 衍生 版 本 . 该 系列 模型 的 
基本 思想 是 将 光 压 摄 动 加 速度 分 解 到 DYB 坐 标 系 
的 3 个 坐标 轴 方 向 ,并 将 每 个 分 量 采用 多 项 式 结合 
周期 函数 形式 进行 描述 . 因此 , 卫星 总 的 光 压 摄 动 
加 速度 asnp 可 表示 为 : 


AU? Fo 

asRP = — — —3 X; (ap: ep + ay: ey + ap: eg) , 
Irs = vr M 

(1) 


式 中 , AU 为 天 文 单位 , AU = 1.495978707 x 10" m; 
rs 和 7 分 别 为 太阳 和 卫星 的 位 置 矢量 ; Fy? 
验光 压 模 型 计算 的 先 验 力 ; M 为 卫星 质量 ; ep, 
ev 和 ea 分 别 为 DYB 坐 标 系 3 个 轴 的 方向 矢量 ， 其 
中 ep 为 卫星 至 太阳 方向 的 单位 矢量 ,ey 为 卫星 至 
地 球 方向 与 卫星 至 太阳 方向 又 乘 的 单位 矢量 , ep = 
ep X ev; ap. ay flag TH] D, YMB 3 个 方向 的 
光 压 模型 分 量 系数 . 

在 BDS 卫 星 精 密 定 轨 中 使 用 的 经 验 型 光 压 
模型 主要 有 ECOMI 模 型 (9 参数 、7 参 数 和 5 参数 )、 
ECOM2 模 型 及 ECOMC 模 型 (13 参 数 )3 个 大 类 . 

ECOM1-9 模 型 ( 即 早期 的 ECOM 模 型 ) 是 目前 
各 分 析 中 心 使 用 最 多 、 定 轨 效 果 较 好 的 光 压 模型 ， 
其 具体 表达 式 为 : 


at 


64 卷 李鹏 杰 等 : 北斗 三 号 卫 
ap = Do + Dc cosu + Dssinu, 
ay = Yo + Yo cosu + Yesinu, (2) 
ag = Bo + Bocosu+ Bgsinu, 

式 中 ,ww 是 卫星 在 轨道 面 内 相对 于 升 交 点 的 角度 ; 


Do. Do. Dg. Yo. Yo. Ys. Bo. Box BsÆ 
ECOMI1 光 压 模型 的 待 估 参 数 . FEA, Do. YoRBo 
4 UD. YMB 3 个 方向 的 常数 项 参数 , Dc 和 Ds 
为 D 方 向 的 一 阶 余弦 周期 项 参数 和 一 阶 正弦 周期 
项 参数 , Yo 和 站 为 Y 方 向 的 一 阶 余弦 周期 项 参数 和 
一 阶 正弦 周期 项 参数 ，Bc 和 Bs 为 B 方 向 的 一 阶 余 
弦 周 期 项 参数 和 一 阶 正弦 周期 项 参数 . 由 于 该 模 
型 9 个 待 佑 光 压 参数 之 间 的 强 相关 性 , 在 实际 BDS 
卫星 定 轨 应 用 时 一 般 采 用 简化 ECOMI1 模 型 ， 即 
ECOM1-7 模 型 和 ECOM1-5 模 型 (51. 

Arnold 等 Ro 为 了 消除 ECOM1 模 型 在 ( 俄 罗 
斯 ) 全 球 卫 星 导 航 系 统 (Global Navigation Satellite 
System，GLONASS) 卫 星 定 轨 中 存在 的 误差 ， 在 
ECOMI1 光 压 模 型 的 基础 上 提出 了 改进 型 ECOM2 
(Empirical CODE Orbit Model 2) 模 型 . 该 模型 具 
体形 式 如 下 : 


ap = Do + ` (Dac cos 2kAu+ 


k=1 
D»ys sin 2kAu) 3 
ay = Yo ; 
QB = Bo + » [Box -c COS (2k = 1) Au+ 
k=1 


Beak—1)s Sin (2k — 1) Au], 


(3) 
XB, Au = u 一 uo0，Uo 为 太阳 在 轨道 面 投 
影 相对 于 升 交 点 的 角度 有 为 正 整 数 ;，np 和 np 


分 别 为 万 方向 、B 方 向 上 的 周期 项 个 数 ; Dakc 和 
Dzxs 分 别 为 D 方 向 的 2k 阶 余弦 周期 项 参数 和 2k 阶 
正 纺 周期 项 参数 ; Btzs_c 和 Bizs_1ys 分 别 为 B 方 向 
的 2k — 1 阶 余 弦 周 期 项 参数 和 2k — 1 阶 正弦 周期 项 
欧洲 定 轨 中 心 (Centre for Orbit Determina- 
tion in Europe, CODE) 通 过 大 量 轨道 数据 测试 , 给 


精密 定 轨 中 的 光 压 模型 


E 
SN 
= 


出 了 ECOM2 (D4B1) 模 型 、ECOM2 (D2B1) 模 型 
两 种 方案 . ECOM2 (D4B1) 模 型 形式 如 下 : 


ap = Do + Dac cos 2A + Dos sin 2Au+ 
Dac cos 4Au + Dag sin 4Au, 


(4) 
ay = Yo , 
ag = Bo + Bc cos Au + Bs sin Au, 
AF, Dac 和 Da2s 为 万 方向 的 二 阶 余弦 周期 项 参 


数 和 二 阶 正 弦 周 期 项 参数 ; Dic 和 D4s 为 万 方向 
的 四 阶 余弦 周期 项 参数 和 四 阶 正弦 周期 项 参数 . 
ECOM2 (D2B1) 模 型 形式 如 下 : 


ab = Do + Doc cos 2A + Dos sin2Aw, 
(5) 


ay — Yo, 


ag = Bo + Bc cos Au + Bs sin Au. 


此 外 麻 省 理工 大 学 (Massachusetts Institute of 
Technology，MIT) 还 基于 大 量 轨 道 数 据 再 分 析 并 
结合 ECOM1 和 ECOM2 模 型 的 优势 , 提出 了 新 的 混 
合 光 压 模型 (Empirical CODE Orbit Model 1 and 
2 Combined, ECOMC), 即 BCOMC 模 型 BI. 目前 ， 
ECOMC 模 型 已 应 用 于 GAMIT 10.7 (GPS Analy- 
sis at MIT, GAMIT) 软 件 和 Ginan 软 件 . 该 模型 的 
具体 表达 式 为 : 

ap = Do + Do cos Au + Dg sin Au+ 
D»c cos 2A + Dog sin 2Au+ 
Dac cos 4Au + Das sin 4Au, 

ay = Yo + Yc cos Au + Ys sin Au, 


(6) 


ag = Bo + Bc cos Au + Bs sin Au, 
日 上 式 可 知 ECOMC 模 型 共有 13 个 光 压 参数 . 


2.1 ” 光 压 模型 的 先 验 力 计 算 

对 于 BDS 卫 星 的 光 压 模型 先 验 力 ( 醒 ) 主 要 通 
过 对 不 同 卫 星 间 面 质 比 系数 进行 换算 、 采 用 简化 
Box-Wing 模 型 计算 等 方式 得 到 . 根据 Box-Wing 模 
型 , 可 将 BDS 卫 星 分 为 若干 个 平面 分 别 进行 处 理 ， 
并 将 各 个 面 的 光 压 摄 动力 进行 合成 , 得 到 BDS 卫 星 
质心 总 的 光 压 摄 动力 22. 该 模型 的 具体 形式 如 下 . 


juni 


64 卷 
+F FE 
为 23-24] 
o 
f; = ZUA; COS 0; | (ai + ôi) ( 
c 
2; COS gn, ; 
对 于 太阳 能 帆 板 ， 
338324). 
o 
f; — Ai cos6; 
C 


2.2 


供 
轨 


于 卫星 星体 ， 平 面 ;受到 的 光 压 摄 动力 户 


Sj 十 zni 


mi) + 
3 


(7) 


PZ BI RO D FE BRN 71 f, 


(a, + 0;) si 十 
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ER 


为 平面 


0 十 p; = 1; 0; 为 平面 ?的 法 向 矢量 mj 与 卫星 至 太 B 
量 s; 之 间 的 夹 角 . BDS-3 IGSO 和 MEO 卫 星 
的 几何 、 光 学 特性 参数 见 表 1， 对 于 该 部 分 内 容 的 


方向 矢 


HY A 


XN 


2 (5 + pi cos.) n,| i 


D 星 在 一 个 天 文 单位 处 的 太阳 辐射 通 
量 、 c 为 光速 、 4 ;为 平面 ;的 表面 H: Qi、 o, o; 2) H] 
i 的 吸收 、 漫 反射 和 镜面 反射 系数 ， 且 Qi + 


见 文 献 [25]. 


表 1 BDS-3 卫 星 的 几何 和 光学 特性 参数 


Table 1 Geometric and optical characteristic parameters of BDS-3 satellite 


Satellite Type PRN? SV Component? Area/m? Absorbed (a)  Difused (6) ^ Refected (p) 
4X3 2.86 0.350 0.000 0.650 
C19, C20, C21, 
-X 2.86 0.615 0.385 0.000 
C22, C23, C24, 
" +Y 3.63 0.135 0.856 0.000 
MEO-CAST C32, C33, C36, 
+Z 2.18 0.920 0.080 0.000 
C37, C41, C42, 
-Z 2.18 0.350 0.000 0.650 
C45, C46 
Solar Panels 10.22 x 2 0.920 0.080 0.000 
C25, C26, C27, +X 1.25 0.200 0.800 0.000 
C28, C29, C30, +Y 3.13 0.200 0.800 0.000 
MEO-SECM? 
C34, C35, C43, +Z 2.59 0.200 0.800 0.000 
C44 Solar Panels 5.40 x 2 0.920 0.080 0.000 
+X 8.496 0.350 0.000 0.650 
+Y 7.560 0.135 0.865 0.000 
IGSO C38, C39, C40 
+Z 4.956 0.870 0.130 0.000 
Solar Panels 17.70 x 2 0.920 0.080 0.000 


1 Pseudorandom noise. 


? Satellite vehicle (including Solar Panels). 


3 X, Y and Z are coordinate axes of aster coordinate system. 
^ The MEO satellite manufactured by China Academy of Space Technology (CAST). 
5 The MEO satellite manufactured by Shanghai Engineering Center for Microsatellites (SECM). 


光 压 模型 的 参数 选择 策略 分 析 


本 文采 用 的 轨道 数据 为 武汉 大 学 分 析 中 心 提 


的 BDS-3 卫 星 精 密 星 历 ， 其 采样 间隔 为 15 min. 
道 数 据 的 时 间 跨 度 为 2019 年 年 积 日 第 1 d 至 2022 


年 年 积 
卫星 光 


佳 组 合 


7EBDS-3 B. BY 


策略 . 


日 第 100 d. 通 
压 参 数 变化 规律 ,以 确定 日 


过 分 析 轨 道 拟 合 获得 的 BDS-3 


8 光 压 模型 参数 最 


L 道 拟 合 中 采用 弧 段 长 度 为 1 d 
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的 卫星 星 历数 据 作为 虚拟 观测 值 进行 动力 学 轨道 
WE, 求解 每 颗 卫 星 的 初始 状态 参数 、ECOMC 模 
型 的 光 压 参数 . 在 轨道 拟 合 中 考虑 太阳 和 月 球 等 
行星 引力 摄 动力 、 太 阳光 压 摄 动力 、 相 对 论 效应 、 
地 球 固体 潮汐 和 海洋 潮汐 等 的 影响 . 轨道 拟 合 中 使 
用 的 地 球 重力 场 模型 为 GOCO06s (Gravity Obser- 
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vation Combination 06s) 模 型 6]. BDS-3 IGSO 和 
MEO 卫 星 采 用 ECOMC 光 压 模型 拟 合 的 光 压 参数 
时 间 序 列 如 图 1 和 图 2 所 示 . 图 中 横 坐 标 为 简化 颂 略 
H (Modified Julian Day, MJD) 表 示 的 时 间 , 纵 坐 标 


单位 为 m . s~?. 


压 模型 先 验 加 速度 (F/M) 的 乘积 ， 
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Fig.1 Time sequence diagram of the ECOMC model parameters for the BDS-3 IGSO satellite 


图 1 为 BDS-3 IGSO ] 


数 项 、 周 期 项 光 压 参数 


HJ RECOMC3E He Be A! af 
的 时 间 序 列 图 和 参数 值 


对 应 的 标准 差 (Standard Deviations, STD). 从 


图 1 中 光 压 参数 变化 规 和 


LE 和 参数 值 的 标准 差 可 以 
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看 出 ， 在 采用 ECOMC 模 型 进行 轨道 拟 合 的 结果 
中 , BDS-3 IGSO 卫 星 的 周期 项 光 压 参数 Dc、Ds、 
Dog. Dacs Das. YoRYs HIE BW fa zE (STD/) 
于 0.004), 且 参 数值 趋 近 于 0, 无 明显 周期 性 变化 趋 
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势 . 因此 , ARED Dss Dass Dacs Dass 用 ECOMC 模 型 进行 BDS-3 IGSO 卫 星 定 轨 应 用 中 ， 


也 和 区 参数 项 对 ECOMC 模 型 整体 的 贡献 相对 较 A ULAR HD. Yo. Bo. Do. Ds. Bo. 
小 . 综合 轨道 拟 合 和 轨道 预报 的 实际 效果 ,在 采 Bs 和 Dzc 8 个 光 压 参数 , 即 ECOMC-8 模 型 . 
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图 2 BDS-3 MEO 卫 星 ECOMC 模 型 参数 时 间 序 列 图 


Fig.2 Time sequence diagram of the ECOMC model parameters for the BDS-3 MEO satellite 
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E 
BE 
= 


图 2 中 给 出 了 BDS-3 MEO-CASTRIBDS-3 M- 报 的 步骤 如 下 : 首先 , 对 给 定 卫 星 初始 状态 下 的 轨 
EO-SECM 卫 星 ECOMC 光 压 模型 常数 项 、 周 期 项 。” 道 参数 进行 积分 , 得 到 卫星 轨道 初始 值 ; 然后 , 将 卫 
光 压 参数 的 时 间 序 列 以 及 光 压 参数 值 的 标准 差 . 旺 精 密 星 历 作 为 虚拟 观测 值 , 建立 卫星 动力 学 轨道 
从 图 2 中 光 压 参数 时 间 序 列 可 以 看 出 ，BDS-3 两 。 ” 拟 合 的 观测 方程 , 并 对 观测 方程 进行 最 小 二 乘 估计 


类 MEO 卫 星光 压 参 数 都 随时 间 呈 现 出 一 定 的 变 得 到 卫星 初始 状态 下 轨道 参数 和 动力 学 参数 的 改 
化 规律 ， 其中， 对 于 光 压 模型 的 党 RR». 正 值 ; 最 后 ， 对 改正 后 的 初始 状态 重新 进行 积分 , 即 
YAB), MEO-CASTE BAD). Bo 参数 波动 可 得 到 卫星 的 拟 合 轨道 和 预报 轨道 

较 大 ， 且 Do 参数 表现 出 了 周期 性 . E 比 MEO- 3.1 ”轨道 拟 合 精度 

CAST 卫 星 , MEO-SECM 卫 星 的 D6。、 Bo 参数 波 动 有 oT Beh eH] A H 
AHN. HERKEN AMIS, BOS Rise ten enit einen aii 
3 两 闫 MEO 卫星 的 Ds、 Dos、 Dis 和 Yo 光 压 参数 。 表 2 中 的 光 压 模型 进行 了 BDS-3 IGSO 和 MEO 卫 星 
值 较为 恒定 (STD 小 于 0.002)， 无 明 吕 周期 性 变化 。 轨道 拟 合 数据 处 理 , 在 数据 处 理 过 程 中 , 使 用 了 相 
趋势 ， 且 参数 值 趋 近 于 0， 因 此 ， 可 以 推断 Ps、 — 同 的 卫星 初始 轨道 参数 以 及 相同 的 其 他 摄 动力 模 
we Dzs、Dss 和 Yc 参数 项 对 ECOMC 模 型 整体 的 页 献 型. 对 于 拟 合 弧 段 长 度 为 24h 的 拟 合 轨道 , 每 5 min 
> R 其 对 轨道 拟 合 结果 产生 的 影响 可 以 忽略 。 输出 一 个 拟 合 结果 . 不 同 光 压 模型 得 到 的 IGSO 和 
Ail. 在 实际 BDS-3 MEO 卫 星 定 轨 应 用 中 ， 可 MEO 卫 星 轨 道 拟 合 结果 精度 见 图 3 和 图 上 


— | y ` MZ AH 5 
= 以 考虑 只 估计 ECOMC 模 型 3 个 方向 上 的 常数 参 从 图 3 和 图 4 中 可 以 看 出 , 使 用 ECOMC 模 型 进 
D XDe WRBOUARERUGHEDO. Ye Bo. Pes — 行 轨道 拟 合 的 效果 最 好 ，BDS-3 IGSO 和 MEO 卫 
Dzc 和 Dac, 即 PCOMC-9 模 型 星 轨 道 拟 合 残 差 在 径 向 (Radial)、 切 向 (Along) 和 
法 向 (Cross) 方 向 上 均 优 于 2 cem. 由 于 ECOM1-5 模 
3 ARNAN 型 和 ECOML7 模 型 缺少 D 方 向 上 的 高 阶 周期 项 ， 


为 了 验证 2.2 节 中 提出 的 BDS-3 IGSO、MEO 使 得 IGSO 和 MEO 卫 星 轨 道 拟 合 结果 的 稳定 性 和 
卫星 ECOMC 模 型 光 压 参数 选择 策略 的 合理 性 , 本 拟 合 精度 相对 较 差 . 相 比 ECOM1 和 ECOM2 模 型 ， 
文选 取 了 2022 年 年 积 日 第 105 d28 28107 d 的 BDS-3 ECOMC 模 型 明显 提升 了 MEO 卫 星 在 径 向 和 法 向 
卫星 精密 星 历 作为 虚拟 观测 值 , 研究 不 同 光 压 模 型 上 的 拟 合 精度 . 由 于 ECOMC 模 型 待 估 参 数 较 多 ， 
AN 对 BDS-3 卫 星 轨 道 拟 合 和 预报 的 影响 , 并 通过 轨道 ” 且 参 数 之 间 存 在 强 相 关 性 , 对 于 IGSO 和 MEO 卫 星 
S ” ” 拟 合 精度 和 预报 精度 两 个 方面 进行 评定 . 由 于 各 天 分别 使 用 本 文 2.2 节 中 ECOMC-8 模 型 、ECOMC- 

的 数据 处 理 结 果 基 本 一 致 , 因此 选取 2022 年 年 积 9 模型 进行 轨道 拟 合 ， 获 得 了 与 使 用 ECOMC 模 型 
FI 28107 d 的 处 理 结果 进行 展示 . 卫星 轨道 拟 合 和 预 相同 的 拟 合 效果 . 


表 2 BDS-3 卫 星 的 光 压 模型 参数 组 合 
Table 2 The combined SRP model parameters of BDS-3 satellite 


Model Name SRP Parameters 
ECOMC Do, Yo, Bo, Dc, Ds, Yo, Ys, Bc, Bs, Dac, Dos, Dac, Das 
ECOMC-8 (IGSO) Do, Yo, Bo, Dc, Ds, Bc, Bs, Doe 
ECOMC-9 (MEO) Do, Yo, Bo, Dc, Ys, Bc, Bs, Dac, Dac 
ECOM1-5 Do, Yo, Bo, Bc, Bs 
ECOM1-7 Do, Yo, Bo, Dc, Ds, Bc, Bs 
ECOM2 (D2B1) Do, Yo, Bo, Bc, Bs, Doc, Des 
ECOM2 (D4B1) Do, Yo, Bo, Bc, Bs, Doc, Des, Dac, Das 
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BDS-3 MEO satellite orbit fitting residuals 
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BDS-3 MEO satellite orbit fitting residuals 
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Fig.3 The statistical results of orbit fitting residuals for BDS-3 MEO satellite 
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g BDS-3 IGSO satellite orbit fitting residuals 下 BDS-3 IGSO satellite orbit fitting residuals 
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图 4 BDS-3 IGSO 卫 星 的 轨道 拟 合 残 差 统 计 结果 


Fig.4 The statistical results of orbit fitting residuals for BDS-3 IGSO satellite 
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轨道 预报 精度 
初始 轨道 和 光 压 模型 参数 选择 对 轨道 预报 影 
响 很 大 ， 而 其 他 的 动力 学 模型 误差 对 轨道 预报 
的 影响 可 以 忽略 不 计 ， 因 此 ， 本 文采 用 中 长 期 
(15 d)、 中 期 (7 d) 和 短期 (1 d) 3 种 模式 进行 BDS-3 
IGSO 卫 星 和 MEO 卫 星 轨 道 预报 ， 以 验证 ECOMC 
模型 光 压 参 数 选择 策略 的 合理 性 . 在 轨道 预报 数据 
处 理 中 , IGSO 卫 星 采 用 了 ECOMC 模 型 和 ECOMC- 
8 光 压 模型 ， MEO 卫星 采 用 了 ECOMC 模 型 和 EE- 


3.2 


ECOMC 
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图 5 BDS-3 MEO] 
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COMC-9 光 压 模 型 , 卫星 初始 轨道 参数 以 及 其 他 摄 
动力 模型 均 相 同 . 预报 轨道 通过 与 对 应 弧 段 独立 
的 卫星 精密 轨道 进行 比较 以 评定 预报 精度 ， 评 定 
指标 为 预报 轨道 与 精密 轨道 互 差 在 径 向 、 切 向 和 
法 向 的 RMS (Root Mean Square) 值 . 表 3 中 给 出 了 
BDS-3 IGSO 卫 星 轨 道 预报 精度 的 统计 结果 . 表 4 中 
给 出 了 BDS-3 MEO 卫 星 轨 道 预报 精度 的 统计 结果 ， 
其 中 各 颗 卫 星 预报 1 d4、7 d 和 15 d 的 轨道 精度 统计 
结果 如 图 5、 图 6 和 图 7 所 示 . 
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'The accuracy of 1 d prediction orbit for BDS-3 MEO satellite 
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Fig.7 The accuracy of 15 d prediction orbit for BDS-3 MEO satellite 


X3 BDS-3 IGSO 卫 星 的 预报 轨道 精度 统计 结果 
Table 3 The statistical results of prediction orbit accuracy of BDS-3 IGSO satellite 


1d prediction/m 7 d prediction/m 15 d prediction/m 


Model Type 
Radial Along Cross Radial Along Cross Radial Along Cross 
C38 2.267 11.612 0.392 16.522 492.203 9.528 38.066 2304.486 140.588 
ECOMC C39 2.038 10.834 0.182 13.132 411.585 11.370 27.472 1953.874 142.261 
C40 2.115 10.855 0.195 14.720 444.769 13.108 32.009 2124.860 151.690 
Mean 2.142 11.107 0.274 14.856 450.735 11.429 32.804 2132.550 144.929 
C38 2.477 12.672 0.363 18.041 532.935 9.709 41.408 2491.409 140.900 
ECOMC-8 C39 1.719 9.120 0.155 11.078 350.222 11.432 23.485 1680.316 142.405 
C40 1.976 10.133 0.187 13.733 417.641 13.000 29.805 1999.992 151.490 
Mean 2.081 10.746 0.252 14.569 440.114 11.459 32.427 2084.110 145.007 


#4 BDS-3 MEO 卫 星 的 预报 轨道 精度 统计 结果 
Table 4 The statistical results of prediction orbit accuracy of BDS-3 MEO satellite 


1d prediction/m 7 d prediction/m 15 d prediction/m 


Model Type 

Radial Along Cross Radial Along Cross Radial Along Cross 

Min 0.034 0.073 0.046 0.622 5.219 8.205 4.246 48.085 99.504 

ECOMC Max 0.089 0.511 0.059 1.749 21.184 7.653 6.424 175.388 94.810 
Mean 0.056 0.262 0.062 1.444 13.085 8.356 5.959 106.105 102.597 

Min 0.033 0.072 0.046 0.626 5.057 . 8.204 4.289 48.32 99.526 

ECOMC-9 Max 0.088 0.505 0.059 1.724 21.424 7.646 6.437 173.149 94.854 
Mean 0.054 0.258 0.062 1.387 13.041 8.353 5.911 105.944 102.594 
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从 表 3 中 的 统计 结果 可 以 看 出 , 对 于 IGSO 卫 星 
中 期 (7 d) 和 中 长 期 (15 d) 预 报 轨 道 , 使 用 ECOMC- 
8 模型 进行 预报 ，7 d 和 15 d 的 预报 轨道 在 切 向 的 
RMS 均 值 分 别 为 440.114 m 和 2084.110 m; 使 用 E- 
COMC 模 型 进行 轨道 预报 ，7 d 和 15 d 的 预报 轨道 
在 切 向 的 RMS 均 值 分 别 为 450.735 m 和 2132.550 m. 
相 比 ECOMC 模 型 ，ECOMC-8 模 型 在 切 向 上 的 预 
报 精度 有 一 定 程 度 的 提升 . 

由 表 4 和 图 5-7 可 得 ， 对 于 MEO 卫 星 短期 (1 d) 
轨道 预报 ， 使 用 ECOMC 和 ECOMC-9 模 型 分 别 进 
行 轨道 预报 ， 两 种 光 压 模型 的 预报 轨道 在 径 向 、 
切 向 和 法 向 的 RMS 均 值 差异 在 1 cm 以 内 ， 对 于 
MEO 卫 星 中 期 (7 d) 和 中 长 期 (15 d) 轨 道 预报 , 使 用 
ECOMC 模 型 进行 预报 ，7 d 和 15 d 的 预报 轨道 在 
径 向 的 RMS 均 值 分 别 为 1.444 m 和 5.959 m、 切 向 
RMS 均 值 分 别 为 13.085 m 和 106.105 m; 使 用 ECO- 
MC-9 模 型 进行 轨道 预报 ，7 d 和 15 d 的 预报 轨道 
在 径 向 的 RMS 均 值 分 别 为 1.387 m 和 5.911 m、 切 
向 的 RMS 均 值 分 别 为 13.041 m 和 105.944 m. FALL 
ECOMC 模 型 ， ECOMC-9 模 型 在 径 向 和 切 向 上 的 
预报 精度 都 有 一 定 程 度 的 提升 . 


4 结论 


本 文通 过 对 BDS 卫 星 精 密 定 轨 中 常用 的 光 压 
模型 进行 研究 , 利用 BDS-3 卫 星 元 数据 ,基于 Box- 
Wing 模 型 计算 了 卫星 的 光 压 模型 先 验 力 . 采用 卫 
星 动力 学 轨道 拟 合 方法 得 到 了 ECOMC 模 型 的 13 个 
光 压 参数 , 并 对 BDS-3 IGSO 卫 星 和 MEO 卫 星 的 光 
压 参 数 进行 了 时 间 序 列 分 析 , 分 别 给 出 了 适用 于 
BDS-3 IGSO 卫 星 和 MEO 卫 星 的 光 压 参数 选择 策 
Wk. 利用 轨道 拟 合 和 轨道 预报 试验 , 对 BDS-3 ISO 


卫星 和 MEO 卫 星光 压 参数 选择 策略 的 合理 性 进行 
了 验证 . 结果 表明 : 
(1) 对 于 BDS-3 MEO LH, 采用 只 包含 常数 参 


数 Do。、 了 Ww 和 Bo 以 及 周期 项 参数 Dc、Ys、Boc、Bs、 
D2c 和 和 Dac 的 ECOMC-9 模 型 进行 轨道 拟 合 , 能 够 获 
得 与 使 用 ECOMC 模 型 相同 的 拟 合 效 果 . 同时 , 使 
jECOMC-9 模 型 进行 中 长 期 (15 d) 和 中 期 (7 d) 轨 
道 预 报 ， 对 MEO 卫 星 轨 道 的 径 向 和 切 向 的 预报 精 
度 有 一 定 的 改善 作用 ; 
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(23K FH RHD Yo. Bo. Do. Dg. Bc. 
Bs 和 Dc 8 个 光 压 参数 的 ECOMC-8 模 型 进行 BDS- 
3IGSO 卫 星 轨 道 拟 合 和 轨道 预报 , 对 于 轨道 拟 合 结 
R, 3 颗 IGSO 卫 星 的 轨道 拟 合 残 差 在 径 向 、 切 向 和 
法 向 上 都 在 2 cm 以 内 . 对 于 7 d 以 上 的 轨道 预报 , 使 
用 ECOMC-8 模 型 的 切 向 预报 精度 有 一 定 程 度 的 提 
FF. 
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A Study of Solar Radiation Pressure Model in BDS-3 Precise 
Orbit Determination 


LI Peng-jie ZHANG Han-wei XIE Meng-xin ZHAO Dong-fang 
(School of Surveying and Land Information Engineering, Henan Polytechnic University, Jiaozuo 454000) 


AnsrRACT For the BDS (BeiDou Navigation Satellite System) satellites, the solar radiation pressure 
is the major non-gravitational perturbation the satellites suffer. Due to the influence of multiple factors, 
accurate modeling the solar radiation pressure acting on these satellites is difficult to accomplish. The 
solar radiation pressure perturbation is an essential source of error in the precise orbit determination 
(POD) and orbit prediction processes of the BDS satellites. As the ECOMC (Empirical CODE Orbit 
Model 1 and 2 Combined) model adopts the advantageous features of the ECOM1 (Empirical CODE 
Orbit Model 1) model and the ECOM2 (Empirical CODE Orbit Model 2) model, a greater number of 
parameters have been introduced into the model, leading to strong correlation between the ECOMC 
model parameters. In response to the deficiency of the ECOMC model, the authors collected the BDS-3 
satellite precise ephemeris data that was provided by Wuhan University (WHU) from January 2019 to 
April 2022. And by these ephemeris data, the ECOMC model’s 13 parameters were obtained through the 
dynamic orbit fitting method. By conducting these time sequence diagrams analysis of the ECOMC model 
parameters, this thesis affords the selection strategies of the solar radiation pressure parameters of the 
ECOMC model for BDS-3 inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO) satellite. 
And the rationality of the selection strategies of the solar radiation pressure parameters for the ECOMC 
model was verified, by orbit fitting and orbit prediction experiments. The results show that using the 
improved ECOMC model gains the best orbit fitting effects for BDS-3 IGSO and MEO satellite and can 
also enhance the accuracy of medium and long-term orbit prediction of BDS-3 IGSO and MEO satellite. 


Key words solar radiation pressure, precise orbit determination, BDS-3, orbit fitting, orbit prediction 
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